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BENDING STRESSES DUE TO TORSION 
IN A TAPERED BOX BEAM 
By Edwin T. Kruszewski 


SLTMMARY 


A method is presented for the calculation of "bending stresses 
due to torsion in a tapered "box "beam. A special taper was assumed 
in which all flanges, if extended, would meet at a point. The general 
procedure of analysis given is similar to the procedure for a non- 
tapered "beam presented "by Raul Kuhn in his paper "A Method of 
Calculating Bending Stresses Due to Torsion, " NACA ARE Dec. 1942. 
Recurrence formulas developed for use in this calculation are 
included- A comparison was made of flange and sheet stresses in 
"boxes with varying taper, including a non tapered "box. 

The results o"btained "by this method were compared with experi- 
mental data o"btained from tests performed on a tapered "box "beam. 

The "box "beam was tested under two independent cond.itions of loading: 
first a concentrated torque at the tip and later, a concentrated 
torque at the quarter point of the span. The experimental results 
o"btained from those tests showed good agreement with the calculated 
results . Calculations for the test specimen for the two loading 
conditions are also shown. 


INTRODLETION 


The basic load -carrying structure of many aircraft ^i^igs is a 
"box of approximately rectangular cross section consisting of 
front and rear spars and the top and bottom skins of the wing. This 
box is stressed by both bending and torsional loads . Only the 
torsional loads are discussed in the present paper. 

The determination of the stress distribution in a box under 
torsion is a relatively simple problem provided that the 
section of the box is not constrained in any manner against w rp g. 
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In this case the well-knovn Bredt formula for thin"walled torsion 
tubes is applicable. If some restraint is offered to varping, 
however, a set of secondary stresses is introduced in the box. 

Because the resuJ tarts of these secondary stresses are actually 
bending moments in the planes of the whlls and are accompanied by 
the shear forces necessary to cause spanwise variation of the 
bending moments, these stresses are usually referred to as bending 
stresses due to torsion. 

A method for the calculation of these bending stresses is 
presented in reference 1 . In order to simplify the calculation, 

Kuhn utilizes an assumption that the cross "sectional dimensions 
and the torqi;es are constant within each bay of the box and gives 
the solution only for boxes in which the sides are parallel. 

Actual wings, however, are usually tapered both in depth and in 
width. The present work is intended to furnish a theoretical 
solution of the effects of taper on bending stresses due to torsion 
and also to present experimental verification of the method of 
calculation. 

The body of this paper is divided into two parts. The first 
part deals entirely vrith the theoretical development of the formulas 
for tending stresses due to torsion in a tapered box beam. The 
general procedure of analysis presented is similar to the procedures 
presented by Kvihn in reference 1 and by-Ebner in reference 2 . In 
order to simplify the mathematical analysis a special taper is 
assumed^ that is, the sides of tho box are assumed to taper linearly 
in such a way as to meet at a point. 

The second part of the paper deals with the experimental 
verification of tho theoretical formulas. A description of the test 
specimen and the test setup is given. Comparisons are then made of 
the calcxilated and experimental results for two independent conditions 
of loading. Complete numerical solutions for both loading conditions 
are given in the appendix. 


SYlfflOLS 

A effective flange area, square inches 

A-r. area of flange angle, square inches 
F 

A area of cover stringer, square inches 

s 
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E Young 's mod^ilus of elasticity, pai 

F flange load at any point, pomde 

G shear modulus, psi 

Kj_, Kp, M taper constants 

L length of flange in individual hay, inches 

/ ^n ^n \ 

E taper ratio [ or 1 

\Vl <=n-l/ 

T external torque, inch -pounds 

V volme of material, cubic inches 

X redundant f 3 .ange force, pounds 

a distance hetveen hulldieads, inches 

h width of cover, inches 

c depth of spar, inches 

f, J, p i^arping constants 

m designation for general term in series 

n designation for typical hay 

q shear flow, pounds per inch 

r designation of root hay 

s distance in each hay measured along axis in flange, inches 

t sheet thickness, inches 

w warping deformation, inches 

X distance in each hay from outboard bulkhead measured 

along axis, inches 

a.. angle between flange and center line of cover, radians 
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angle tetveen flange and center line of spar, radians 
a normal flange stress, psi 

T shear stress, psi 

average value of shear stress, psi 

Suh scripts; 


h refers to covers 

c refers to spars 

i desiignates inboard end of hay 

o designates outboard end of bay 


Superscripts: 


T designates stresses di;e to torque 

U designates stresses due to dummy unit loads 

X designates stresses due to X -forces 

The subscripts of the redundant flange force X and of 
dimensions b and c designate the stations at which they exist, 
whereas the subscripts for T, w, p, J, and f designate the bay 
under consideration. 

The bulkheads or stations are denoted by 0, 1, 2,... n-1, n, 
n+1, . . . r, starting from the tip or outboard end and proceeding to 
the root or inboard end. (See fig. 1.) The bays are also numbered 
from the tip, the tip bay being designated as number one (see 
fig. l) • A bay therefore carries the nimiber of its inboeird bulkhead 
or station. 
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IJEVELOPMEKT OF TBEOKETICAL FOEMULAS 


In actual wing design neither the spanwise vaxiaticn of the 
torque nor the cross “sectional dimensions can he represented hy 
simple mathematical expressions . In order to simplify the 
mathematical calciilations the box is divided at the bulkheads into 
a n’jmber of bays end the torque is assumed constant within each 
individual bay. 


A box beam under torsion is an indeterminate structure. In 
order to make the structvire statically determinate, the box is cut 
at each bulkhead and redundant flange forces X are applied at 
all flanges. These X -forces are axial forces applied at each 
flange, as shown in figure 2. Under the action of the torque T 
and the flange forces X the box is deformed as shown by the 
dashed lines in figure 2. The amount of deformation is calculated 
by the use of the principle, of virttial work, sometimes knovni as 
the dummy -unit -load method. The X -forces are then found by the 
application of the principle of consistent deformation of adjacent 
bays . 


Siggi conve n tio n. - External torques T are positive when acting 
clockwise as viewed from the tip. Shear stresses t are positive 
when acting in the direction of shear stresses caused by positive 
torque. The X-forces are positive when acting in the directions 
shown in the sketch in figure 2. Normal stresses a are positive 
when caused by positive X-forces. The warping deformation w is 
positive in the direction as shovm by the dashed lines in figure 2. 

General assuiautions . - The cross section of the box is assumed to 
be rectangular and doubly sjTEmetrical. The shape of the cross section 
is maintained by the bulkheads, which ai-e assumed to be rigid in 
their own planes. In place of the actual structure, the equivalent 
structiuce shown in figure 2 is used, in which all the area capable 
of carrying normal stresses is concentrated in the flanges. The 
walls of the equivalent structiure are assumed to carry cnly shear 
stresses and the flanges, all the normal stresses . In order to 
allow for the fact that the walls can actually carry normal stresses, 
each flange area is increased by one -sixth of the area of both a 
cover and a spar web. If the covers include stringers, an effective 
stringer area is added to each flange. This effective stringer area 
is that area which, when concentrated at the flange, gives the same 
section modulus about the neutral axis of the cover as the actual 
stringers. In the case of equally spaced stringers the effective 
stringer area is simply one -sixth of the total stringer area. The 
taper of the box is such that the flanges meet at a point. The 
effective flange area is assumed constant within each bay. 
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Stresses in an Individual Bay 


The formulas given herein are derived for a typical hay n, 
hounded hy the hulhheads n-1 and n, (fig. 2). The hay is acted 
on hy three independent sets of loads : a torque T^^ on both ends 

of the hay, a group of -forces on the inhoard end, and a group 

of X , -forces on the outboard end. Formulas are derived for the 
n~l 

stresses due to each of these independent loads. The final stress 
distribution may he obtained hy superposing’ the individual stresses . 

Stresses caused hy toroue . - Tlie shear stresses caused hy the 
torque acting on a hay are given hy the well-}cnovm formula for 
shells in torsion (Brodt's formilLa) 


■where h is the v/idth of the cover and c the dep'th of the spar 

at some distance x from the outboard end (fig* 3) * Iii 'the case of 

pure torque nc normal sti'esses are set up in the flange. 

Str esses caused hy the X~force s.~ Vvhen a set of X -forces is 
applied~to the end of the hay, both axial stresses in the flanges 

and shear stresses in the vails are set up within the hay. Unlike 

the shear stresses for the nentapered tox, the shear stresses in tho 
tapered box are not constant throughout the hay . 

In order to study the distribution of the shear stresses 
within the individual wells, a section of a wall is isolated as 
shown in figure 3(a). The free tody sho-wn in this figure is a part 
of a spar web hounded by two planes cutting the spar just inside 
the flanges and hy two parallel planes, one just inside the^ n-3. bulk- 
head and another parallel, to it and at a distance x from it. The 


1 


T 


h ^ 2hct^ 


( 1 ) 


T 
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loading on the tody is also shown in the figure . By simmation of 
moments about the point of intersection of the flanges, an expression 
is obtained for the shear flow in the spar in terms of the out- 
board shear flow , 

°n-l 


<1 


c 



Q r» , 

“'^n-1 


(2a) 


and, similarly. In the cover 



(2b) 


where the notation is the same as that shovm in figure 3(^)* 

The distribution of the shear flow q„ is obtained from an 

infinitesimal section of a spar web isolated as ^ho'wn in figure 3(b). 
The free body sho'vm in this figure is a section of the spar web 
bounded by two parallel planes an ijrfinitesimal distance dx apart 
and by two plaiies cutting the spar just inside the flange. From 
a summation of moments, the fundamental shear-flow relation in the 
spar for a tapered box beam is obtained. 


'Ic 



( 3 a) 


and, similarly, in the cover 


% = % 

X 

This equation shows that at every point along the box the shear flow 
acting on the flange is equal to the shear flow in tlie walls at 
that point. 

In order to obtain a relationship between the shear forces in 
the cover and in the spars, a free body of a cross section of the 
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"box is considered. Since the taper of the "box is such that the 
center lines of the flanges meet at a point, the flange loads 
contrihute no torque end the condition 2T = 0 gives the equation, 


q-jjhc + q^jbc = 0 


or 




i^) 


Two expressions for the f3.ange loads are derived, one for the 

-forces and another for the X -forces. The free body of the 
n n-1 

flange in figure 3(c) shows the loading of the flange under X^ -forces 

only. A summation of forces along the flange sho\na in figure 3(c) 
gives an expression for the flange load F at any point along the 
bay: 


F 


ds 

dx 


;'x 

i 

I 

Jo 



ds , 


nx 


■0 


q-b dx = 0 

X 


(5) 


where — — is a constant depending on the taper of the box. The 
dx 

combination of equations (2) to ( 5 ) gives 


F 




ds 

dx 



( 6 ) 


If the integration is performed with c given as a function 
of X by the equation 


"n-1 


1 




3C 

a 


c 


(7) 
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equation (6) tecoiaes 


F = 2L q 


Ct,- 


n-1 X 


®n-l c 


a 


( 8 ) 


where L is the length of the flange and a is the distance hetween 
"bulkheads . 

VTith the forces X applied at the outboard end of the hay, 
the flange loads hecoime 


^ = ^n-1 + 'ic 1- 

°n-l c a 


( 9 ) 


In order to obtain q in terms of the X -forces, the value 

^n-1 

of F for both X . - and X -forces is calculated for x = a. 
n-1 n 

For the case of the Xj^ -forces 


X = 2q 

""n-l Cn 


or 


X. 


Ic " "'1‘h 
°n-l 


= E 


n 


h-1 


2L 


( 10 ) 


For the case of the X -forces 

n-1 


^n -1 

0 = X T + 2q„ L 

n-1 - 


"n-1 Cn 


or 


X. 


n-1 


= 


= -B 


n-1 


n-1 

2L 


( 11 ) 
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t>n 

where R is the taper ratio = 

l^n-1 


n 


°n-l 


By substituting equation (10 ) in equations { 2 ), (4), and (8) and 

equation (11) in equations (2), (U), and (9), a summary of tho stresses 

in the bay in terms of the X -forces at the end can be made. For 

the X -forces 
n 


= °n °n -l in 
r.2 2L 


■') 


S 




X 


n “n-l __n 

,2 2l 


> 


(12) 


0 



c a A 


and for the X , -forces 
n-1 


<1 


°n^n-l ^n-1 

o2 21, 





(13) 


where A is the effective area of the flange. 
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Deformation of an Individioal Bay 


Princ l nle of calcu].atlon . - Under the action of the torque 
and groups of X -forces the cross section of the box warps out of 
its plans, as shc'.im in figure 2. The magnitude of this warping is 
calculated by the method of virtual work. The following three steps 
are necessarj'’ to obtain the magnitude of the warping. First, the 
stresses T and a due to the ^plied loads arc obtained and 
second, the stresses and c due to a system of dummj'’ loads 

are calculated. These dumiTiy loads are unit loads applied at the 
point where the deflections are desired and also in the direction 
desired. The last step is to obtain the defomation by use of the 
equation 


where V is the volume of the stressed material. This equation is 
based on the principle that the external work done by the unit force 
must equal the internal energy stored by virtue of the existence of 
tho unit force. 

Examination of figure 2 shows that the warping is doubly 
antisymmetrical end consequently the duamy loads emp3.oyed in the 
solution are also doubly anti syumetri cal. This group of unit loads 
is similar to the group of X -forces and therefore the formulas for 
the X -forces can be used in the calculation of the stresses caused 
by dumi^y loads . 

W arning caused by torque .- The stresses caused by the torque 
acting on the box are, from formula (l), 



1 





( 15 ) 
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In order to obtain the warping in the nth hay at huUdiead n, the 
anti symmetrical group of unit loads is applied at the inhoard end 
of the hay. The stresses caused hy these forces are calculated 
from fo3rmulas ( 12 ) hy placing = - 1 . The results are 


U 

a = 


^n X 1 
c a A 


U 


c c 1 1 

n n -1 

• .2 2 Lt. 







(16) 


' T ^ = Vn-1 1 

21 , j 

The results for a, and and for 0^^, """c^ 

given in equations (1?) and (16) are now substituted in equation ( 14 ) 
to give 


4 w, 


n,- 


= 2 


Jo ^ 2 hct(^ 


Cn‘^n"l 


c 


2 



dXc 


+ 2 



2: ^n I ^n^n ~1 

G 2hctb \ t)2 



htb dx^ 


which yields when integrated 


V 

" 8G\.-l^n-l \ ^ H 


] 

L % J 


1 + R 


2 E 


(17) 
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In order to ottaJ.n the varping due to torsion at hulkhead n-1, 
the groups of unit loads are applied at the outboard end of the bay. 
Now the stresses due to the dummy loads are calculated from 
formula (I3) with = !• The results are as follows: 





> (18) 


I 

j 


The stresses due to the torque acting on the bay are those 
shown in equation (15)- An inspection of the stresses (equations (16) 
and (18)) caused by the unit loads gives 




(19) 


Warping caused by X^-forces.- The warping at the inboard end 
of the bay caused by the -forces is obtained by ‘the application 
of equation (l4) in the form 



U 

= Ij. / A ds + 2 

Jo E 


Jo 


btjj dXb 


+ 2 


/■■a. T T U 

/ cTqT(, 


Jc 


G 


ct dx 
c c 


( 20 ) 


Substitution of the values for the stresses due to the unit loads 
(equation (16)) and the stresses due to the X,^-forces (equation (12)) 
in equation (20), causes the equation for warping to become 


X. 


n 


w. 


^i 


= -Vn 


( 21 ) 


Ik 
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where 


Pn 





Pn-1 + ^ °n-l 
t-b L tc 


^ + B 
2 



and 




(B - 1)3 


(^1 . a 10s, K 


( 23 ) 


The constant is dependent entire 1;^ on the taper ratio for 

the individual hays. 

As the va].ue of the taper ratio approaches unity, the 
expression for in equation ( 23 ) hecoraes too sensitive for 

practical use. By expanding the i.ogerithm into an infinite series 
of (B - 1), 

K, = 1 + - (B - I) - - (B - 1)^ + i(B - 1)^ .... 

-‘■2 5 10 


+ 


m+1 

("D 


35 

(m + l)(m + 2)(m + 3) 


m 

(R - 1) 


(24) 


This expression can only he used for small values of B, since the 
series converges onlj’’ for values of B < 2. The rate of convergence 
is very slow, however, after B reaches a value of approximately 1.5* 
For values of B < 1.5 only four terms of the series are needed to 
evaluate K 3 within 1 percent. A graph of numerical values for 

from B = 1.00 to B = 3.00 can he seen in figure 4. 


The warping at the outboard end of the hay caused hy the 
Xjj “forces can he written in the form 



( 25 ) 


= -Jn^n 
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where the coefficient is ohtained by the application of 

equation (l4) . By substituting the values for unit stresses from 
equation (18) and the valv.e for X^-s tresses from equation (12) into 
equation (l4) and integrating the result, is found to be given by 


6ae 


1 / Q-b ^n-T ^ ^ ^ n-1 \ B + 1 

8gL \ L t^) L j 2 


(26) 


where 



2E^ - 2R - (1 + E) E logg E 
(E - 1)3 


(27) 


For small values of E, again the expanded series for 


^2 


= 1 - i (E > 1)^ + ^ (E - 1)^ . . . 


10 


10 


+ (- 1 ) 


m+1 6(m - 1) 


m 


(m + 1) (m + 2) (m + 3) 


(E - 1) (28) 


is more practical. As for K^, the rate of convergence of the 

series makes the expression practical only for the range from E = 1.00 
to E = 1.50. The numerical values for are plotted in figure 4. 


Warping caused by -forces . - The warping at the inboard 
end of the bay, caused by the X^_^ -forces, can be shown, by means 
of Maxwell’s law of reciprocal deflection, to be equal to 


^n-1 


w 

^i 


"^n^n-l 


(29) 


where is the expression given in eqxiation (26) . 
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The '.Tarping at the oufboerd end of the hay is derived hy 
Btihstitutlng the values of the stresses for imit Idads from 
equation (18) and the values of the stresses for the -forces 

from equation (13) into equation (l4). Upon integrating the result, 
the follo-^ing expression is obtained; 


Xn-1 ■ 

^^n^ = Vn-1 


(30) 


where 


f ^ ^n-1 ] 1 + B 

n 3 3 EA 8gi, y L t^ L tc 


(33.) 


and 


K, = 3 


- 1) - 2R log E 


(32) 


For small values of E, the expanded series form for K^, 


Kj = 1 - i- (B - 1) (E - if (E - 1)3 . . 


*(■ 1 ) 


m 


(m 4- 2)(m + 3) 


(P - 1) 


m 


(33) 


should te used. As vas the case for the expressions for and Kg, 

the rate of convergence of the series inahes the expression for 

practical only within the limits of E = 1»00 and R = 1.50» 

Comparison of tanered-hox formulas with m iform. cross "section 
formu las , - A comparison of the formulas derived for the tapered 
i)ox~heams with those formulas derived for the non tapered hox 
hy Kuhn in reference 1 can easily he made hy obtaining the formulas 
for the special case of taper where the taper ratio is unity. 
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\Jith this assumption, the following relations hold true 


a = a^ = = L 


Ki = Kj = Kj = 1 


Equation (I7) now hecomes 


■ 8Gbc 


T T 



which is identical with equation (21) of reference 1. Equations (22) 
and (31) become 


which is identical with equation (23) of reference 1. Equation (26) 
becomes 


which is identical with equations (27) of reference 1. 

This comparison shows that all the tapered “box formulas for the 
special case of E = 1.0 or uniform cross section are identical 
with those formulas developed in reference 1. 

Derivation of a recurrence formula . - The total warping at the 
ends of the bay n due to the combined effects of the three separate 
forces T, and Xn-1 can now be obtained by a svamnation of all 

the component parts . The equation for the warping at the inboard end of 
bay n is given by the sum of equations (I9), (21), and (29) 




T 


(3i^) 


^n^^ - ■'^n " Pn^n '' '^nr^n-1 
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The equation for warping at the outboard end of bay n given by 
the euraniation of equations (I9), (25), and (30) is 


The equation for warping at the outboard end of bay n + 1 is 
obtained merely by increasing all subscripts of equation ( 35 ) "by 
one and is 


T 

'^(n+l)Q " '^n+l " <^n+l^n+l ^n+l^n 


(36) 


According to the principle of consistent deformation, the 
warping at the outboard end of one bay must be equal to the warping 
at the inboard end of the adjacent bey. A recurrence formula can 
therefore be obtained by equatiiig the warping formulas of two 
adjacent bays. By equating the expression for w^^ in equation ( 34 ) 


and the expression for 


^(n+1), 


in equation ( 3 ^c), the general 


recurrence formula becomes 


Jn\-l ■ I ^n+1 * % j ='n + ‘ "n* (3T) 


By giving n successive values from n=l to n=r, a 
set of r equations is obtained, each of which contains three of 
the redundant X -forces. These equations represent the continuity 
conditions at stations 1 to r. The tip of the bo:: is usually free 
from any restraint, thereby. making the force Xq at the tip of the 

box zero. Therefore, for a box divided into r bays, r rodundant 
forces and r equations exist. 

•Boundary condit ion. - With the tip- of the box free from any 
restraints and the tip force X^ = 0, the first equation of the 

system is 


- ^ f2 + Pi I Xi + ~ ^^2^ 



(38) 
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VHien a team is attached to a rigid foundation, the fovindation may 
te considered a tay r + 1 having infinite shear stiffness and 
infinite axial stiffness; therefore. 


^V+3. ~ ^r+1 "^r+l “ ^ 

The last equation of the system now tecomes 

Vr-3. ■ Vr =■ V (39) 


: . Comnarison of stx'*esaes for tox e s with varying' taner ratio and 
flange area . - In order to show, the effect of taper on the tending 
stresses due to torsion, calcula.tions were made of tlie flange and 
sheet stresses for toxes with varying taper ratio and flange area. 

All toxes were 120 inches long and w'ere divided into six hays of 
equal length. The root sections of the toxes were identical. The 
dimensions for the toxes at the root aaid tip can te seen in tatle 1. 

In figure 5(a), curves of the flange stresses ere dra'vm for 
the three boxes with constant flange areas and also for the 
moderately and highly tapered toxes with a tapered flange area. 

In figure 5(^)^ curves of the sheet stresses are drawn only for the 
three toxes with the ta.pered flange area. 

Figure 5(s-) shows that the flange stresses for the case of the 
moderately tapered box (curves C and D) are only slightly greater 
than those of the nontapered box (ciurve E) • Examination of the root 
flange sti’esses for all cases shows that the increase in these 
stresses for an incx’ease in taper is very small; the root flange 
stress for the highly tapered box is approximately 10 percent above 
that of the nontapered box. Since, as in the nontepered box, the 
flange stresses for the box with modLerate taper decrease very rapidly, 
the only appreciable flange stresses occur in the vicinity of the root. 
As the taper tecomes .greater however, the flange stresses along the 
total span increase . Inspection of curves A and B shows that 
the flange sti’esses at a point 20 inches from the tip are approxl " 
mately one -third of the maximum flange stresses at the root. 
Consequently, for the box team with moderate taper the tending stresses 
due to torsion in the outboard end are negligible as compared with 
the stresses near the root, whereas in the highly tapered box the 
flange stresses do become appreciable. A similar conclusion for the 
sheet stresses can be reached by the inspection of figure 5(t)* 


20 


HACA TW Wo. 1297 


Comparison of curves A and E and curves C and D in 
figvire 5(a) shows the effect of a varJatjon of flange area along 
the span on the flange stresses . The maximum effect appears near 
the tip of the hox where the flange stress for the box with tapered- 
flange areas is approximate 20 percent higher than the flange 
stress in the box A^ith constant flange area. 

Two sets of calculations for the root stresses were made for 
a moderately tapered box under a distributed torq,ue loading. The 
distribution of the torque was such that the increment of torque 
in each bay was proportional to the chord of the bay at the inboard 
end. The method presented herein is used for the first calculation 
which was an exact soli’.tion of the root stresses. The second 
calculation was made on the assumption that the box consisted only of 
the root bay with the total torque acting on the outboard end. The 
result of these caD.culations showed that the root stresses calculated 
by means of the exact m^ethod were approximately 10 percent greater 
than the root stresses calculated by the approximate method. Similar 
results A^ere observed for a calculation of the stresses in the highly 
tapered box under a similar distribution of torque. VJhen only the 
approximate value of the stresses at the root is required, a satis- 
factory answer can be obtained by assuming that the box consists only 
of the root bay and that the total torque is concentrated at the 
outboard end of that bay . 

In both sets of calculations A.'ith the distributed torque, the 
flange stresses near the tip of all the boxes \rere approximately of 
the same order of magnitu.de as those of the moderately tapered box 
with tip loading shoAm in figAjre 5 (s ) • 


EXEEP.IMENT/.L i’ERIFICATION OF FOPi'iUIi^S 


Test snecimens . - In order to obtain experim.ental verification 
of the formulas derived herein, a large tapered box beam was 
constructed and tested. The box A^as made sj/mmetrical about the midspan 
and was supported there by a rigid frame as shown in figure 6. 

Because of this setup, complete restraint against A-rarping at the 
root, which is the midspan of the box, was ass-ured. Equal torques 
were’ applied at the tips and later at the quarter points of the span. 

The material used for the box was 24s-T aluminum alloy and the 
general dimensions of the specimen are shoAra in figure 7. The 
thicknesses of the covers end the spar A/ebs and the sizes of the 
flanges, stiffeners, and stringers are also shoAm in figure 7 * 

General dimensions and stringer spacing of Doth the root and 


MCA Oil No. 1297 


21 


tip sections are shown in figure 8, along with the dimensions for 
the equivalent substitute sections. The dimensions of the box 
vary lineai’ly from the tip to the root, whereas the thicknesses 
of the sheets and sizes of the flanges and stringers are constant 
for the whole box. Although the flange area of the simplified 
structure varies linearly from the tip to the root, due to the 
addition of one -sixth the area of the cover and spar web, an 
average flange area in each bay is used in the calculations . In- 
order to assTire that the bulkheads were rigid for all practical 


purposes, the bulkheads were made of formed ^-inch steel plate. 


Test setun . - The general test setup is shorn in figure 6. 

The box beam was connected to the center tov/er by means of foxir 
steel flexure plates, one on each side of the box. In order to 
reduce the end effects as much as possible, the box was connected 

to the flexure plates at the root by closely spaced ^-inch bolts. 

This type of connection permits the torque reaction to be distributed 
as a uniform shear flow around the perimeter of the box. Fig\ire 6 
shows the loading arrangment at the tips of the box. The same 
method of loading was used when the tests were run with loads at 
the quarter -span points. 


Test procedure . - Strain siurveys were made for both loadings 
with 2 -inch Tuckerman optical strain gages. Shear-strain measurements 
wore taken around the perimeter of the box at sections near the 

center line of each bay, and also at sections l| inches on either 

side of bulkheads 4 and 5 . The shear-strain measurements across 
any cross section consisted in measurements made between the 
stringers and between the flange and adjacent stringer on the covers 
and three equally spaced measurements in the spar web. The shear 
stresses were obtained from strain readings at ky and 135 "tli® 

axis of the structure. The normal strains were measured along each 

flange at approximately intervals starting at 4^ inches from 

the root. 


For each test run, strain -gage readings were taken at zero 
torque and after each of four equal Increments of 75,000 inch-poiinds 
of torque . The load ws.s then released and another zero reading was 
taken. If the two sets of zero readings did not agree within 100 psi, 
a new test run was made. The strain readings for each gage were then 
plotted against torque and the best straight line was drawn through 
the points. If the line did not intersect the origin, a parallel line 
was drawn through the origin. If, however, the new lino was displaced 
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from the original line "by more than a strain equivalent to 2C0 psi, 
a new set of readings was taken. The values of strain were then 
obtained from this new parallel line. 

Test results . - In converting the results obtained from the 
strain surveys to stresses. Young's modulus was taken as 10,6o0.ksi 
and the shear modulus as 4,000 ksi. 

The observed shear stress of a cover or a spar at any section 
is the average of the shear stresses obtained for the two opposite 
covers or spars at that section. Figure 9 shows tlaat the distribution 
of the shear stresses across the covers and spars is constant 
througliout the section, except for the part of the skin between 
the flange and the adjacent stringer. These plots substantiate in 
pajrt the original assumption of uniform shear stresses over the 
cross section. 

The observed flange stress at any section is the average of 
the flange stresses obtained for the four flanges at that section. 

In figvire 10 a plot of the stringer stresses at a cross section 
in bays 5 s.nd 6 is shovn. Strain readings were taken on each 
stringer on the leg adjacent to the cover. The cross section in 
bays 5 ^d 6 were chosen because the normal stresses were the largest 
in those bays. Figure 10 shows that the chordwise distribution of 
the stringer stresses in the cover is approximately linear; thereby 
the use of the theoretical equivalent -area coefficient of one -sixth 
appears to be justified. 

Comparison of test results with theoretical curves . - A comparison 
of the observed and calculated shear flow and normal flange stresses 
for both loading conditions can be seen in figures 11 to l4. The 
stresses were calculated as shown in the appendix by means of the 
formulas presented in this paper. 

Examination of figures 11 and 3.2 shows very good agreement 
between the observed and calculated shear -stress values. The only 
point in the tip loading condition (fig. 11) that does not fall 
on the calcxilated curve within the accuracy of the Tuckerman gage 
readings is a point on the spar 7-5 inches from the tip, which is 
approximately 5 percent greater than the cozrresponding calculated 
value . This deviation from the calculated ciorve is probably caused 
by the fact that the section at vrhich this measurement is taken is 
near the tip of the box where the load is being applied. 

Examination of figiu’es I3 and.l4 shows good agreement between 
the experimental and calculated values of the axial loads in the 
flange due to torsion. For both loading conditions, the observed 
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values in lays 5 ''■^hich include the most critical loads (the 

loads at the root of the hox), agree within the accur’ecy of the 
Tuckerman gage reading. The observed values in bays 3 compare 

favorably with the calculated values for the tip -loading condition. 
For the quarter -point loading condition, however, the observed 
values for the test points in the vicinity of bullchead 3 a^re 
slightJ.y greater than the calcxxlated values for the seme points • 

This deviation from the calculated values can be explained by the 
uncertainty of the conditions at the point of loading. The extent 
to which the loading fixture restrains the flanges from warping 
and the fact that the torque applied to the box through the spar 
webs needs an appreciable distance to be distributed, are only the 
obvious reasons for deviations to occur in the vicinity of the 
loading buJ-kheads . Also for these reasons, no attempt was made 
to evaluate the stresses in the tip bay for the tip loading condition. 

The values of the experimental stresses in figiares 13 and l 4 , 
taken within 5 inches of bulkheads 1 and 2 , are somewhat greater 
than the calculated values . This deviation may be explained in 
part by local bending of the flange. The flange acts as a 
continuous beam supported by an elastic foundation and loaded at 
each bulkhead. The bending moment in this beam is the greatest 
at the bulkheads where the deviation of the observed values from the 
calculated values are the greatest. 


CONCLUDING REMARKS 


The agreement that was obtained between the experimental 
bending stresses due to torsion in a tapered box beejn and the 
calculated values indicates that the method presented can be used 
to obtain the bending-stress distribution in a tapered box beam 
imder torsion. 

For boxes with very small taper the flange stresses in the 
outboard half of the box are very small. In these boxes a first 
approximation of the bending stresses due to torsion can be made 
by using the properties of the tapered box at the root, by 
considering the box nontapered, and then using the method of 
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analysis descrl'bed "by Paul .Kuhn In his paper on tending stresses 
due to torsion (NACA AER, Dec. 1942). For mere accurate calculations 
of the tending stresses due to torsion in a tapered tox, however, the 
method presented herein shoiAld te used. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., February 13, 19^7 
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APIENDIX 

SOIUTIOIJ 07 STRESS DTSTPI3TTTI0N IN TAPERED BOX BEAM 


In order to give an illustration of the method, a complete 
solution of the shear and normal stress distribution in the 
tapered "box used as a test specimen is given. 

The actual cver-all dimensionc of the specimen are given in 
figure 7. The dimensions at the root and tip cross sections are 
given in figure 8, together witlx the dimensions for the equivalent 
structure- The effective flange area is obtained from th.e 
eqixation 


A =A^+4-^A„ + ^ bt. + ct^ (Al) 

F o S 6 D o 2 


vhore A_, is the area of a flange angle and Ac- is the area 

of covGi* striiigor* Tho effective flange areas are assumed "to 
be concentrated at the points of intersection of the center lines 
of the cover sheets and spar webs. All properties of the equivalent 
box taper linearly from tip to root. Tho box is divided into 
six bays. The geometrical properties for the Inboard and outboard 
end of each bay are listed in table 2. 

The values of E and G used in these calculations are, 
respectively, 10,600 ksi and U,000 ksi. For the first loading 
condition, a torque of 100,000 Inch-pounds is applied to the tip 
of the box and for the second loading condition a torque of 
100,000 inch-pounds is applied to the middle bulkhead of the box. 

The warping constants p, and f, calculated by means 
of the equations (22), (26), and ( 3 I), a-^d the waiu^ing due to 
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torsion calculated ty the use of equation (1?); tabulated 

in the following table. The tabulated w^ values are for the load 
at tip of box. 


Bay 

Warping constants 

Pn 


fn 

T 

1 

1.610 X 10 

-6 

0.829 X 10 

/• 

1.564 X 10 ^ 

-6 

6439 X 10 

2 

1.581 

■ 574 

1.511 

5848 

3 

1.647 • 

.684 

1.585 

5274 

4 

1.714 

.794 

1.660 

4803 

5 

1-785 

.904 

1-739 

4410 

6 

1.857 

1 1.011 

1.817 

4083 


For a box with six bays, the recurrence formula (equation ( 37 )) 
and the equations for the boundary conditions (equations (38) and ( 39 )) 
give the following six equations. 


" ( ^1 ^2)^1 <) 2^2 “ ^ 2 *^ " 

JgX3_ - (p2 + f3]x2 + J3X3 = W3'^ - W2^ 

'^ 3^'2 “ (^3 )h " ^" 3 ^ > (A 2 ) 

J 4 X 3 - (P4 + f5 )X4 + J5X^ = w^T . ,,,T I 

j5Xj^ - (P5 + f6 )^5 ^6X5 = ■ ^'^ 5 ^ I 

J6^5 ■ P6 ^6 = "'^ 5 ^ 
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From substitution of the values from the preceding table into the 
formulas of equation (A2), the following set of simultaneous equations 
is obtained for the first tip loading condition; 


^1 

X2 

^3 

^4 

='5 

"^6 

■W . T - 

n+1 n 

-3-121 

.574 

0.574 
-3 • 166 
.684 

0.68i+ 

-3.307 

.794 

c .794 

-3.453 

.904 

0.904 

-3-602 

1.011 

1.011 

-1.857 

-591 

-574 

-471 

-393 

-327 

“1 k )83 


The solution of these equations gives 


= 2 hk Xj^ = 437 

Xg =292 x^ = 965 

x^ = 307 Xg = 2724 

These values give the flange loads at each bulkhead. In order to 
obtain the distribution of the flange loads between the bulkheads, 
the formula 




(A3) 


obtained from a summation of equations (12) and (I3) is used. 

For calculation of the distribution of the shear flow along 
the span of the box, the formulas 


n^ ^n-l _ 1 - X T ^ ^ 

c2 2L ' ^ 


> (Ak) 


X 




T 

2b c 


and 
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ottained from a summation of equations (1), (12), and (13) are used. 

A tabulation of the flange loads and shear flows can be seen in 
table 3 • 

For the quarter -point loading condition the values for the 
warping constants p, J, and f are the same as those for the 
tip loading condition. Also the values for w^^^ are the same except 

in bays 1 , 2 , and 3 where w^^ is zero. The formulas (A 2 ) can still 

be used for this loading condition. By substituting the appropriate 
values in equation (A 2 ), the following set of simultaneous equations 
is obtained for the quarter-point loading condition: 


1 — 1 

X 

^2 1 

^3 

X4 

"5 

^6 

'’^n+l " '"n 

-3.121 

0.57*t 

-3.166 

.684 

0.684 

-3.307 

.79^^ 

1 

0.79^ 

-3.t^53 

.904 

0.904 

-3.602 

i.oii 

1.011 

-1.857 

0 

0 

4803 

-393 

-327 

"!i083 


The solution of these equations gives 

= -63 = -21 

Xg = -3h3 X 5 = 829 

X = -1528 = 2650 

3 ° 

Again, by the use of the equations (A3) and (a 4 ) the distribution, 
of the flange loads and shear flows was obtained along the span. 
The tabulations of these flange loads and shear flows are shown 
in table U. 
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TABI£ 1 


EOOT AMD TIP DIMENSIC®S OF BOX BEA^B 


Section 

b 

(in.) 

(iS.) 

(in.) 

tc 

(in.) 

A 

(Ccaistant) 
(sq in. ) 

A 

(Tapered) 
(eq in.) 

Root 

50.000 

10 .000 

o.o 64 

0.072 

1.000 

1.000 

Tip: 

Highly 

tapered 

20.000 

4.000 

.o 64 

♦ 

.072 

1.000 

.200 

Moderately 

tapered 

35-000 

7-000 

.064 

.072 

1.000 

.200 

Non tapered 

50 .000 

10.000 

.o 64 

.072 

1.000 

.200 


NATIONAL ADVISORY 
COMMITTEE FOR AERONAUTICS 


8 


NACA TN No. 1297 


TABIi 2 


CO 


GEOMETRICAL PROEEBTIES OF BAYS IN TAPERED BOX BEAM 


Bay 

(1^.) 

®'C 

(In.) 

% 

(in.) 

L 

(in.) 

‘^n-l 

(In.) 

°n 

(in.) 

^n-l 

(In.) 

^n 

(In.) 

tc 

(in. ) 

tb 

(in.) 

A 

(eq in.) 

(0^) 

B 

^1 


^3 

1 

15,0 

15.050 

15.000 

15 .050 

5.869 

6 . 4 l 6 

26.132 

28.602 

0.0700 

0.0629 

0.935 

1.093 

1.045 

0.998 

0.955 

2 

20.0 

20.067 

20.003 

20 .071 

6 . hl 6 

7.146 

28.602 

31.896 

.0700 

.0629 

.972 

1.114 

1.054 

•997 

.946 

3 

20.0 

20.067 

20.003 

20.071 

7.146 

7-^5 

31.896 

35.190 

.0700 

.0629 

1.015 

1.102 

1.050 

.998 

.950 


20.0 

20.067 

20.003 

20.071 

7-^5 

8.605 

35.190 

38.484 

.0700 

.0629 

1.059 

1.093 

1.045 

.998 

.955 

5 

20.0 

20.067 

20.003 

20 .071 

8.605 

9 . 33 ^^ 

38.484 

41.778 

.0700 

.0629 

1.101 

1.085 

1.040 

•998 

.960 

6 

20.0 

20.067 

20.003 

20 .071 

9.331^ 

10.064 

41.778 

45.072 

e0700 

.0629 

1.145 

1.078 

1.037 

.999 

.963 


*ATorei6® of offcctlre areas at Inboard and outboard ends of bay. 
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TABLE 3 


CALCULATED NOEMAL AND SHEAE STRESSES IN 
TAIEKED BOX BEAM FOE LOAD AT TIP 
[t = 100,000 in. -It] 


Bay 

Dlstaince 
from tip 
(in.) 

( 11 >) 

0 

(psi) 

% 

(It /in.) 

% 

(psi) 

% 

(It/in.) 

(psi) 


0 

0 

0 

33^^.9 

4780 

317.1 

504 o 

1 

10 

167 

178 

296.7 


280.9 

4470 


15 

2 lf 4 

257 

279.9 

4ooo 

265.1 

4210 


15 

244 

257 

273.8 

3910 

271.2 

4310 

2 

25 

269 

277 

245.0 

3500 

242.6 

3860 


35 

292 

294 

220.4 

3150 

218.2 

31^70 


35 

292 

294 

219.7 

311*0 

218.9 

31^80 

3 

lf 5 

300 

296 

198.9 

2840 

198.1 

3150 


55 

307 

296 

180.7 

2580 

180.1 

2860 


55 

307 

296 

183.9 

2630 

176.9 

2810 

If 

65 

375 

35 ^ 

167.9 

2400 

161.5 

2570 


75 

^37 

405 

154.0 

2200 

148.0 

2350 


75 

h37 

405 

165.3 

2360 

136.7 

2170 

5 

85 

712 

647 

152.0 

2170 

125.8 

2000 


95 

965 

859 

140.3 

2000 

116.1 

1850 


95 

965 

859 

175-3 

2500 

81.1 

1290 

6 

105 

1878 

l64o 

162.3 

2320 

75.1 

1190 


115 

2724 

2336 

150.8 

2150 

69.6 

1110 
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TABLE It 


CALCULATED NOEMAL AND SHEAE STBESSES IN TAEEBED BOX 
BEAM POE LOAD AT QUAETEE POINT OF SPAN 


Bay 

Distance 
from tip 
(in.) 

F 

(lb) 

a 

(pel) 

(Ib/in.) 

7 'c 

(l>si) 

1 b 

(Ib/in.) 

(pel) 

T = 0 


0 

0 

0 

-2.3 

-30 

2.3 

40 

1 

10 

-lt 3 

-46 

-2.0 

-30 

2.0 

30 


15 

-63 

-66 

-1.9 

-30 

1-9 

30 


15 

-63 

-66 

-7.8 

-no 

7.8 

120 

2 

25 

-211 

-217 

-6.9 

-100 

6.9 

110 


35 

- 3'^3 

-345 

-6.3 

-90 

6.3 

100 


35 

- 3^3 

-345 

-32.5 

- 46 o 

32.5 

520 

3 

lt 5 

-965 

-951 

-29.4 

-420 

29.4 

470 


55 

-1528 

-1473 

-26.8 

-380 

26.8 

430 




T = 100,000 ln.“ 

It 




55 

-1528 

-1473 

221.5 

3160 

139-3 

2210 

It 

65 

-740 

-699 

202.1 

2890 

127-3 

2020 


75 

-21 

-19 

185.4 

2650 

116.6 

1850 


75 

-21 

-19 

174.0 

2490 

128.0 

2030 

5 

85 

421 

•382 

160.0 

2290 

117.8 

1870 


95 

829 

738 

147.8 

2110 

108.6 

1730 


95 

829 

738 

177-0 

2530 

79-4 

1260 

6 

105 

1774 

1549 

163.9 

2340 

73-5 

1170 


115 

2650 

2273 

152.3 

2180 

68.1 

1080 
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Figure I. — Convention for numbering 
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(a) Free-body sketch of section of spar web. 




c 




c+dc 




dx 
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(b) Infinitesimal section of spar web. 


(c) Free-body sketch of flange. 


Figure 3.— Free-body sketches of component ports of 


'^par. 
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Figure 4.— Variation of taper constants with taper ratio 
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Fig. 4 NACA TN No. 1297 





Flange stress, psi 


NACA TN No. 1297 


Fig. 5a 



Distance from tip of box, in. 
(a) Flange stresses. 
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Figure 5 — Comparison of stresses in boxes with varying toper. 
T = 100,000 inch-pounds. 


Additional shear flow due to 
flange loads (qc~ T/2bc), in-lb 



(b) Sheet stresses. 
Figure 5. — Concluded. 


Fig. 5b NACA TN No. 1297 
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Fig. 6 



Figure 6.- Setup for torsion tests of tapered box beam. 
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Figure 7.— Tapered box beam. 
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(a) Actual root and tip cross sections. 



Root Tip 
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(b) Equivalent root ond tip cross sections. committee aeronautics 


Figure 8. Actual and equivolent root and tip cross sections. 
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Figure 9 — Distribution of shear stresses over cross section 


T= 100,000 inch-pounds. 


CD 
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Section B 
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Figure 10 — Distribution of stringer stresses over cross section. T-lOO, 000 inch pounds. 


Fig. 10 • NACA TN No. 1297 


Shear flow. Ib/in. 
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Fig. 11 
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Figure II. — Comparison of cxporimcntol and calculated sheor flows for tip looding condition. 

T 'lOOpOO inch-pounds 
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Shear flow, Ib/in. 


Fig. 12 
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Calculated Bredt shear flow 

Calculated spar shear flow 

Calculated cover shear flow 

n Experimental spar shear flow 



t)istance from tip of box, in. 
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Figure 12 — Comparison of experimental and calculated shear flows for quarter— point 
loading condition. T= 100,000 inch-pounds. 
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Fig. 13 
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Figure 13.— Comparison of experimental and calculated flange loads for tip loading condition. 

T= 100,000 inch-pounds. 


Flange load, lb 


Fig. 14 
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Figure 14. — Comparison ol experimental and calculated flange loads for quarter- point 
loading condition. T= 100,000 inch-pounds. 


